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Recent experimental observations have found two different kinds of "strange kinetic behaviors" 
in individual semiconductor nanocrystals (or quantum dots). Fluorescence intermittency observed 
in the quantum dots shows power-law statistics in both on and off times. Spectral diffusion of the 
quantum dots is also described by power-law statistics in the sojourn times. Motivated by these 
experimental observations we consider two different but related problems: (a) a stochastic lineshape 
theory for the Kubo- Anderson oscillator whose frequency modulation follows power-law statistics 
and (b) photon counting statistics of quantum dots whose intensity fluctuation is characterized by 
power-law kinetics. In the first problem, we derive an analytical expression for the lineshape formula 
and find rich type of behaviors when compared with the standard theory. For example, new type 
of resonances and narrowing behavior have been found. We show that the lineshape is extremely 
sensitive to the way the system is prepared at time t — and discuss the problem of stationarity. In 
the second problem, we use semiclassical photon counting statistics to characterize the fluctuation of 
the photon counts emitted from quantum dots. We show that the photon counting statistics problem 
can be mapped onto a Levy walk process. We find unusually large fluctuations in the photon counts 
that have not been encountered previously. In particular, we show that Q may increase in time even 
in the long time limit. 



I. INTRODUCTION 



Recent studies have found fluorescence intermittency 
phenomena in single semiconductor nanocrystals (or 
quantum dots) such as CdSc illuminated under a con- 
tinuous wave laser field [Q-Q. In these experiments, a 
quantum dot (QD) typically exhibits blinking behavior; 
at random times the QD jumps between a bright state in 
which it emits many photons and a dark state in which 
it is "turned off" Q-|7]. The on (or off state) is believed 
to correspond to a single electron-hole pair (or ionized) 
state of the QD Thus, the statistics of on and off 
times can tell us about the kinetic mechanisms of the 
QD blinking process. 

In dramatic contrast to the usual expectation, distri- 
butions of on and off times of QDs follow a universal 
power-law behavior, not the characteristic, exponential 
behavior of Poissonian kinetics. It was found that the 
probability density functions (PDFs) of on and off times 
decay as P on (t n) oc l/C on and P oS {t oS ) oc 1/Q ft , 
where rn on ( ff) « 3/2 [Q-wj. The off time distributions 
were measured over more than five decades in time, and 
seven decades in the PDF, -P ff(^off)- PJ-@|. This behav- 
ior appears universal; it is found in all individual QDs 
investigated, independent of the temperature, radius of 
the QD, and the laser intensity J§-£|. The on time dis- 
tributions exhibit similar features. Although a secondary 
photo-induced mechanism introduces a cut-off time in the 
PDF of on times, the power-law behavior has still been 
observed in on time statistics over four decades in time 
and five decades in the PDF, P on {t n) |§-§|- The above 
mentioned cut-off time in the on time PDF depends on 
the laser intensity and the temperature, and when these 
effects become small, the cut-off time appears to diverge 



(i. e., power-law behavior with no cut-off time). A single 
computational realization of the intensity fluctuation of 
the QD based on the two-state model is shown in Fig. [l]. 

The spectral diffusion process of the QD has also been 
investigated in other studies In these studies, the 

fluorescence emission spectrum of QDs was found to fluc- 
tuate between two central frequencies, for example, u+ 
and oj-. The statistics of times, t+ and i_, during which 
a quantum dot emits photons at the frequency of u> + 
or lu^ was studied, and it was found that PDFs for t + 
and t_ follow the power-law behavior, P±{t±) oc 
Moreover, the exponents m± are similar to those in the 
blinking statistics, m± w 3/2 p0[ . This suggests that 
there exists a strong correlation between fluorescence in- 
termittency and spectral diffusion in QDs ||[l0|]. 

Standard approaches to lineshape phenomena and pho- 
ton counting statistics are based on the Markovian as- 
sumption, and the spectral fluctuation is characterized 
by a finite, microscopic timescale. The statistics of t on 
and t g (or t + or £_) observed in QDs clearly indicate 
the breakdown of standard assumptions made in conven- 
tional lineshape and photon counting statistics theories. 
The dynamical behaviors of the spectral diffusion pro- 
cess and the blinking process in QDs is non-Mar kovian, 
non-stationary, and non-ergodic. All these features are 
related to the fact that the average on and off times di- 
verge, for example, (i ff) oc J °° dtt~ 3 / 2 t = oo. Hence, 
the dynamics of spectral diffusion and intermittency ki- 
netics in QDs cannot be characterized by a microscopic 
timescale, which results in "strange" kinetics and leads 
to "strange" conclusions. 

The i -3 / 2 power-law behavior indicates that a simple 
random walk mechanism may be responsible for the ob- 
served behavior. Consider the following scenario: The 
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electron-hole pair generated in the QD under illumina- 
tion is ionized via various mechanisms such as thermal 
H or Auger ionization and the ionized electron or 
hole performs a certain kind of one-dimensional random 
walk either in a physical space such as on the surface of 
the QD, or in the energy space. Then, as is well known, 
the PDF of the first return time of the random walker to 
the origin follows the mentioned 3/2 power-law behav- 
ior This PDF corresponds to the off time PDF. 
Similarly, the on time PDF will exhibit a 3/2 power- 
law behavior when the ionization event is controlled by a 
one dimensional random walk mechanism. Shimizu et al. 
have suggested a resonant tunneling mechanism ||. In 
this mechanism, the energy level of the ionized electron- 
hole pair in the QD in the dark state performs a one- 
dimensional random walk process, and when it matches a 
certain resonant condition, the tunneling process is facil- 
itated so that it enables the recombination of the ionized 
electron-hole pair to occur, thus "turning on" the QD 
Alternatively, Kuno et al. have proposed a mechanism 
that attributes the power-law intermittency to fluctua- 
tions in the environment where QDs are located Q . Lo- 
cal changes in the QD environments may cause the width 
and height of the tunneling barriers for the recombina- 
tion of the electron-hole pair to fluctuate. If the Wentzel- 
Kramers-Brillouin-type theory for the tunneling process 
is considered, a minor fluctuation in the width or height 
of the tunneling barrier may result in a broad range of 
the tunneling times j^]. However, this does not explain 
the observed, universal 3/2 power-law behavior. As far 
as we are aware, there is no definite physical picture of 
the exact nature of these processes. 

Motivated by these observations, we consider in this 
paper two different but closely related phenomena char- 
acterized by the same power-law stochastic processes: 
lineshape theory with spectral diffusion process and pho- 
ton counting statistics of fluorescence intermittency. As 
we demonstrate below, the lineshape for an ensemble 
of systems with a power-law spectral diffusion process 
such as QDs is completely different than that with the 
usual Poissonian case. An important issue in a power- 
law stochastic process is stationarity, which is of concern 
due to a very broad temporal distribution for underlying 
processes. We take into account the stationarity issue in 
the lineshape problem, and show that the lineshape is a 
very sensitive measure of stationarity when the underly- 
ing dynamics obeys power-law statistics. 

In the second problem, we consider the photon count- 
ing statistics of QDs undergoing fluorescence intermit- 
tency characterized by the power-law process, and show 
that it exhibits unusually large fluctuations of photon 
number counts not encountered previously. We will show 
the relation between the statistics of photon emitted from 
QDs and Levy walk processes that have been introduced 
in the context of continuous time random walks p3| . [l4[ . 
Levy walk processes have been used to describe an en- 
hanced diffusion (i.e., super-diffusion), and applied to 
many cases including a tracer diffusion in rotating flows 



Pq] , models of deterministic chaotic diffusion [[L6| and 
of diffusion in random environments fl7| , [l8|| . We pre- 
dict that the photon counting statistics of an ensemble 
of these systems will exhibit large deviations from ordi- 
nary photon counting statistics. Specifically, Mandel's Q 
parameter that measures the fluctuation of the photon 
counts will increase as the measurement time increases 
even in the long time limit. 



II. LINESHAPE THEORY 

Since its introduction by Kubo and Anderson (KA) 
|]l9| , f20f in the context of the lineshape theory, stochastic 
approaches to spectral lineshape theory have found wide 
applications in condensed phase spectroscopy ranging 
from magnetic resonance spectroscopy [19 20 , nonlin- 
ear spectroscopy pl|-p4| to single molecul 
^5r 



c spectroscopy 

Analysis of lineshapes observed from an en- 
semble of molecules as well as from single molecules 
have revealed important dynamical information on the 
interaction between the chromophore and the environ- 
ment. When molecules are embedded in a condensed 
media, the absorption frequency of the molecules changes 
in time due to the interaction between the molecules 
and the environment, which leads to a spectral diffusion 
process [ ^lpfj| -|33| . The influence of the spectral diffu- 
sion on the lineshape has been studied in many cases 
|||l 25 27],|34|,|35|] , and one of the well-known examples in 
these studies is the motional narrowing phenomenon: the 
linewidth decreases as the bath fluctuation rate increases. 
The motional narrowing phenomenon has been observed 
at the level of both the ensemble t j36pl| ] and the single 
molecule |37|| . In this section, we consider the Kubo- 
Anderson oscillator whose frequency undergoes spectral 
diffusion characterized by the power-law process as ob- 
served in recent studies of QDs |],[l(| , and calculate the 
lineshape of the oscillator. 



A. Kubo-Anderson Oscillator Model 

The stochastic lineshape theory of Kubo and Ander- 
son is based on the equation of motion for the transition 
dipole (T9]|2(]|||, 



dt 



(2.1) 



where tu(t) is the stochastic frequency of the oscillator. 
The dynamical quantity which determines the lineshape 
is the dipole correlation function or the relaxation func- 
tion defined by 



*(t,t ) = (n(t)n(t )*), 



(2.2) 



where the average is taken over all the possible realiza- 
tions of the underlying stochastic process and we have 
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set (Im^o)! 2 } = 1- From the equation of motion we can 
calculate the relaxation function as 



(2.3) 



$(*,to) = (exp [i / drw(r) 



When the process is assumed to be stationary, 
<J>(t,to) =<&{t — to), then the normalized steady-state line- 
shape I(u>) can be calculated as the Fourier transform 
of the relaxation function by making use of the Wiener- 
Khintchine(WK) theorem p§], 



i r°° l 

= — / dte-^it) = -Re$(iw + e), (2.4) 

27T J_„ 7T 



where the symmetry of $(t), i) = $*(£), has been 
used, and e — > + . The Laplace transform of z(t) from t 
to s has been denoted by 



z(s) = C{z{t)} 



Ate- St z(t). 



(2.5) 



To derive the lineshape formula given in Eq. (2.4) one 
can start, for example, from the optical Bloch equa- 
tion governed by a stochastic Hamiltonian coupled to 
a monochromatic laser field and use a standard per- 
turbation approximation for the field-matter interaction 

In calculating the average in Eq. ( |2.3| ) , we assume that 
the underlying stochastic process is a renewal process 
as in the KA approach. To make the model as simple 
as possible, we only consider a two state model in this 
work. However, the current formulation can be extended 
into multi-state case. The transition frequency w(t) of 
the chromophore can take the value of either uj + or w_ 
depending on the perturber state, + or — , respectively. 

Each alternating path between the states + and — of 
the perturber leads to a stochastic realization of chro- 
mophore frequency modulation, and it is characterized 
by a sequence of sojourn times in the states + and 
— . The sojourn times in the states ±, t±, are as- 
sumed as mutually independent, identically distributed 
random variables described by the probability density 
functions (PDFs), ip±(t±). With the Markovian assump- 
tions the original KA process amounts to the exponential 
sojourn time PDF, 



. . . 1 f t± 
V±(t±) = — exp 



(2.6) 



We do not assume any specific functional forms for the 
sojourn time PDFs from the beginning, but are mainly 
interested in the process where the sojourn times are dis- 
tributed with long time power-law tails, t± 1+a \a > 0). 
Recent studies on the statistics of the spectral diffusion 
process in QDs correspond to the case a ~ 0.5 [ M . 



B. Calculation of Lineshape 

Unlike the Markovian, KA process, care should be 
taken in order to consider the stationarity of the non- 
Markovian stochastic processes described by the so- 
journ time PDF with non-exponential forms. When the 
stochastic process has been going on for long times be- 
fore the beginning of the measurement at time t = 0, it is 
legitimate to assume that stationarity has been achieved 
in the process if (t) is finite. 

We introduce f±(i), PDFs for times at which the os- 
cillator makes the transition ± — * =p for the first time 
after the beginning of the measurement, knowing that 
the oscillator was at ± at t = 0, and these PDFs might 
be taken different from ip±(t) in general. Following the 
physical argument given by Feller and others p^-445|l (see 
also Appendix), the sojourn time PDFs for the first tran- 
sition event after the measurement beginning at t = are 
given by f±(t±) for the stationary case, 



/±(*±) = — /"dr^M 

T± Jt + 



where the mean sojourn times r± are given by 



dtti/j±(t), 



(2.7) 



(2.8) 



and they are assumed to be finite. This type of the 
initial condition is called the equilibrium initial condi- 
tion [43] ], and they should be used in describing the 
stationary stochastic process. In this case, we have 
$(t, to) — Q(t — to), and we set to = 0. When the mean 
sojourn time diverges, the concept of stationarity breaks 
down. 



For the Poissonian case given in Eq. (2.6), we have 
f± = 4'±- Therefore, stationarity is naturally satisfied 
in the Poissonian process. However, the non-Poissonian 
process in which we are interested will not be stationary if 
we simply set f± = ijj±, and the WK theorem therefore 
does not hold in general. We note in passing that sig- 
nificant difference between stationary and nonstationary 
cases has manifested itself in many physical problems, 
e. g. transport properties in disordered materials K3 ] 
and power-spectra in chaotic systems [ [f6| . 

The conditional relaxation functions $ mn (m, n = 
+, — ) are defined over the stochastic paths that start 
from the state m at time and end with the state n at 
time t, 

Ui) = (exp (i [ d™(r) N )\ , (2.9) 

\ \ JO J I mn 

and the total relaxation function is given in terms of & m n 



$(t) 



E 

m— ± n— ± 



Pm$rnn{t)- 



(2.10) 



A sketch of the calculation of Q mn is given by using 
the convolution theorem of the Laplace transform [pTj. 
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For simplicity, we only consider in details the stochastic 
paths of the perturber that begin at the state + at time 
and end at the state — at time t, thus contributing 
to Along a particular path if no transition is ever 

made until t, then the contribution of this path to $++ is 
given by F + (t + )e luJ+t+ in the time domain, where F±(t±) 
the arc the probabilities that the first events ± — * T do 
not happen until time t and given by 



with the initial distribution of the perturber state given 
by 



p± = 



(2.20) 



F±(t ± ) = / dr/±(r). 
Jt± 



Then from Eq. (p]j) the lineshape is given by 



I(w)= -Re 

7T 



(2.11) 



\{-- 









l iV(i-W(i-i) 



This contribution will amount to F + (s — ilu + ) in the 
Laplace domain. The next possible paths are those which 
make the first transition to the state — at t\, jump back 
to the state + after remaining at the state — for time ti , 
and stay at the state + until time t. The contribution of 
these to $ ++ (f) is given by 



T + + T_ \ Z + Z_ 



1 - -0+^- 



, (2.21) 



where z± = ilo — iu>± and ip± — ij)±{z±), and we have 
expressed all the sojourn time PDFs and survival proba- 
bilities in terms of , 



dt 1 



dt? 



dt 3 / + (t 1 ) e i;w +* 1 V-(i2)e^-' 2 * + (t 3 )e^ +t3 (2.12) 



with the constraint t% + t-i + t-j = t. Here, *&±(t±) are 
the survival probabilities corresponding to ip±(t±), and 
defined by 



*±(s±) = 
f±{s±) = 
F±{s±) = 



1 ^ Vj±(f±) 
s± 

1 

T± 

l-/±(*±) 



(2.22) 
(2.23) 
(2.24) 



*±(*±) = / dr^±(r) = r±/±(t±). 
Jt± 



(2.13) 



In the Laplace domain this contribution will read as 

/+(s - ilj + )V3_(s - iw_)* + (s - iu + ) (2.14) 

by the convolution theorem. Summing all the possible 
stochastic paths, we have 

$++(*) = F + (s+) + f + (s+)j>-(s-) 

x [i + v- + ( s+ )v-_( s _) + (^ + ( s+ )^_( s _)) 2 + ■ • ■] 

x *+(s+) 



= F + (s+) + 

i-iM»+)iM«-) 

where s± — s — iuj±. In a similar way, we have 

; + ( S+ )#_( S _) 



$+_( S ) 



(a) = F_(s_) 
$_+(,s) - 



i-V>+(s+W-(0 

/_( S _)V-+( S+ )*_( S _) 



/_( S _)* + ( S+ ) 



(2.15) 

(2.16) 
(2.17) 
(2.18) 



Eq. ( 2.21 ) is the final expression of the lineshape function 
for the stochastic oscillator undergoing the stationary two 
state frequency modulation. 

It is our aim here to show that the lineshape theory ex- 
hibits a very strong sensitivity on the choice of PDF for 
the first event. This becomes important for experimental 
situations when it is not always clear if the underlying 
process is stationary or not. For this purpose we de- 
fine a quasi-lineshape I(u>) by considering the following 
situation: In some experimental situations, a stochas- 
tic process undergoing in the system is not an on-going, 
stationary process, but it is initiated at a certain time, 
for example, t = 0. For this nonstationary case, we as- 
sume that the same sojourn time PDFs, ip±, describe the 
statistics of all the sojourn times, regardless of the first 
or the next sojourn times after t = 0. Then, a relaxation 
function can be defined as Eq. ( |2~3"| ), $(t) = $(t, 0), 
with f± replaced by ip±. The quasi-lineshape I{oS) is 
mathematically defined as a com plex Laplace transform 
of l>(i) in the same way as Eq. (|2.4|). It is obtained by 



replacing f± by ip± in the derivation of Eq. (2.21), which 
yields, 



I{uj) = -Re 

7T 



-P- 



1 



i-V+ , iMi-V»-) 



ip+ip- { \ z + 
-i>- , i>-0-- i>+) 



The total relaxation function can be calculated from 
the conditional relaxation functions, 



(2.25) 



$00 = X! X! P™®"in(i 



(2.19) 



m— ± n— ± 



Note that for Poissonian case I(u) = /(w). However, as 
we show here, a strong sensitivity on the first event is 
exhibited for power-law processes such that I(oj) ^ 
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We emphasize that I{ui) is not a lineshape obtained via 
the WK theorem. In general, the relaxation function, 
$(ii,t 2 ), will depend on two-times for a nonstationary 
process, t\ and t 2l as given in Eq. (2.3), and the corre- 
sponding lineshape will be given by 

J(w,T)~i [ T dt, [ T d^e-^-^^h), (2.26) 
1 Jo Jo 



which will depend on the total measurement ti me T . 
When the stationarity condition is satisfied, Eq. ( 2.26 ) 
is reduced to Eq. (2.4) in the limit T — > oo via the WK 
theorem |B8|. 



C. Examples and Discussion 



The original KA model is recovered from E q. (2.21 ) by 
choosing an exponential sojourn time PDF Jl9|, 20| . We 
first consider the sojourn time PDFs which have finite 
first moments, r± < oo, but divergent second moments. 
As a representative of this class, we use the following 
form, 



1/2 



03/2(i;r)=(^ 



(2.27) 
(2.28) 



In this case, ip±(t±) decays as t ± 5 ^ 2 at long times, thus 
the first moment exists, but the second moment diverges. 

In Fi g. 3 we have compared I(u>) and I(lu) for ip±(t) 
in Eq. ( |2.28| ) and for the KA case. For simplicity, we set 
the average frequency between lu+ and as zero which 
amounts to a simple shift of the frequency origin, and 
define the magnitude of the frequency modulation as uiq, 



LJ + + 10- 



0. 



UJ 



LO + — LU- 



(2.29) 
(2.30) 



We have chosen rp + (t) = ip-(t) = ip 3 / 2 {t; r) in Eq. ( |2.28| ) 
by setting t + = t_ = r. The correlation time of the 
perturber dynamics is varied from slow (luqt 3> 1) to 
fast (luqt <C 1) modulation cases in (a)-(d). For the KA 
case, the well known phenomenon of motional narrowing 
is shown: in the slow modulation case we see two peaks at 
w = ±luq while in the fast modulation case we observe a 
single peak at u> = 0. For the case of i/'3/2) the stationary 
and nonstationary cases show very different behaviors as 
the correlation time is decreased; thus, the first event 
in the underlying random process has a strong effect on 
the lineshapes. In addition, new phenomena are found 
for the stationary lineshape in the fast modulation cases 
(Fig. ^(c) and (d)): three distinct peaks are observed for 
the stationary case. 

The new peaks we observe in Fig. ||(c) and (d) at 
lu = ±luq for the stationary lineshape result from the 



first event in the stochastic process u>(t). The proba- 
bility for the perturber remaining at the initial state is 
governed by the long time tail in the sojou rn t ime PDF. 
Due to the stationarity condition in Eq. (2.7) the sur- 
vival probability for the first event decays more slowly for 
the stationary case(~ i 1_Q ) than for the nonstationary 
case(~ t~ a ), where a — 3/2 in this example. Therefore, 
the stationary case effectively requires the perturber to 
remain at the initial state until much longer times than 
the nonstationary case, resulting in the enhanced peaks 
at uj = ±cjJq. This is why we observe new peaks not 
present in the standard Poissonian case. 

As the next example, we consider the one-sided Levy 
density as the sojourn time PDF Q, 



Tp(s) = L a (rs) = exp(-(rs) Q ), 



(2.31) 



with < a < 1, and r being a coefficient with a time 
dimension. It is well known that the Levy PDF in 



Eq. (2.31) decays algebraically at long times t/r 1, 
L a (t/r) ~ t _ ( 1 + a ) j and thus all the moments of L a (t/r) 
including the first moment diverg e Q. Therefore, there 
is no microscopic timescale for this PDF, and the form of 
f±(t) given in Eq. (2.7) cannot be applied. However, in 
realistic situations, the power-law statistics is modified 
at long times to various reasons, for example, lifetime of 
a molecule. Therefore, it is natural to introduce a cut-off 
time t c such that the algebraic decay is valid during time 
interval r ct <i c . We introduce an exponential cut-off 
function for the convenience of an analytical treatment. 
Now the sojourn time PDFs are given by 



^±(t)=A4e- t / tc L Q (*/r ± ), 



(2.32) 



where Af± are the proper normalization constants de- 
pending on the cut-off time. Then the Laplace domain 
expressions of ip± (t) can be written as 



ip±(s) = exp 



r f) fi - a + <r} 



(2.33) 



Then f±(t) are given from Eq. (2.7) with the mean given 
by 



r± = at c 



(2.34) 



Note that in the limit t c — > oo the Levy PDF without 
cut-off is recovered as ip±(s) = exp(— (r±s) a ) and t± di- 
verge. 

In Fig. |^ we have investigated the effect of the cut- 
off time in the Levy PDF case both in the stationary 
and the nonstationary cases. When a = 0.3 [Figs. | (a) 
and (b)] , both the stationary and the nonstationary Levy 
lineshapes show distinct peaks at u) = ±ujo- As the cut- 
off time is increased, lineshapes become narrower. When 
a = 0.8 [Figs. H (c) and (d)], there appears a peak near 
lu = not present in a = 0.3 case in Figs. |^ (a) and 
(b) in addition to two resonance peaks. This is a new 
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type of the narrowing behavior in that it is controlled by 
the power-law index a rather than the correlation time r 
(as in the Poissonian case), and is termed power-law nar- 
rowing behavior. Also, as the cut-off time is increased, 
the central peak in the lineshape for the stationary case 
diminishes while it remains in the nonstationary case. 
This is because in the stationary the cut-off time 

is increased the first event will dominate the probability 
weight in the stochastic paths of the perturber dynamics. 
The difference between the stationary and nonstationary 
cases is therefore more significant in the Levy case than 
in the case 1/13/2 (i) given in Fig. 1. 

To investigate the power-law narrowing behavior we 
consider the limit t c — > 00. In this limit the stationary 
lineshape approaches two delta functions, 



cos(7ra c /2) 



a c = 0.594611 



(2.41) 



I{tjj) = p+S(u> + lo ) + p_6(w - w ), 



(2.35) 



since the second term in the Eq. ( [2.21 ) vanishes as 
t± — > 00. The nonstationary case, however, yields in 
the limit \cu\r± <C 1, 



We also note that in a recent study of the statistics of 
persistent events that models spin flips separated by ran- 
dom time intervals described by the Levy law it has been 
shown that the distribution of the mean ma gneti zation is 
described by an expression similar to Eq. 
symmetric case, that is, 77 = 1 in Eq. ( 2.40 ) 




for the 



To confirm this finding we have plotted the stationary 
and nonstationary lineshapes for the Levy PDF in Fig. |] 
for the symmetric case (77 = 1) as a is changed. Since 
we have chosen a very large but finite value of t c (= 10 4 ) 
the stationary lineshape still shows the central peak de- 
pending on the value of a, although it will approach two 
delta functions as t c — > 00. The nonstationary case in 
Fig. |^ (b) shows the concave-to-convex transition at the 
critical value of a as predicted. For the stationary line- 
shape, similar kind of behavior can be observed in Fig. ^ 
(a), however, a c now depends on t c . 



lim I(lu) — 

C sin(Tra) 2+x+x- 1 I I ^ 

) 2-kuiq r)x" + (rjx a )- 1 +2 cos(ira) ^ ^ ^° 

I \oj\ > ujq 



(2.36) 



This expression has been obtained from Eq. (2.25) by tak- 
ing the small frequency limit of the sojourn time PDF, 



ij)±(s) = 1 -A±s a 







(2.37) 



where A± = r±. Note that Eq. ( 2.36 ) is not limited 
to Levy PDF, but valid for any PDF with <-( 1+Q ) tail 
(0 < a < 1). Here, a dimensionlcss frequency, x, is de- 
fined by 



LOO +Ll> 

X = , 

LUq — UJ 

and an asymmetry parameter, 77, by 



17 



lim 

t c ^oo p 



p+ 



lim — 

t c — >oo T_ 



A. 



(2.38) 



(2.39) 



Eq. (2.36) shows very asymmetric power-law singulari- 
ties, I(oS) ~ l/(wo i lu) 1 ^" when |cj| < ujo depending on 
a. It is worthwhile to mention that such a strong asym- 
metric lineshape has been encountered in the problem of 
the X-ray edge absorption of metals KTfl. 

In the symmetric case (77 = 1) Eq. ( |2.36|) reduces to a 
simpler expression, 



lim I(lu) = 



sin(7ra) 



2-kloq x a + x 01 + 2 cos(7ra) 



(2.40) 



and zero when \u>\ > ujq. In this case there exists a criti- 
cal value of a below which the lineshape is concave and 
above which convex at to — 0, which is given by 



III. PHOTON COUNTING STATISTICS 

Photon counting statistics has proved useful for inves- 
tigating dynamical processes of an ensemble of molecules 
as well as of single molecules in condensed phases 
pfl| , p8|j29[ . In a semiclassical theory of photon count- 
ing statistics, when there is no source of fluctuation in 
the dynamics of chromophores other than shot noise due 
to the discrete nature of photons, the counting statis- 
tics of the photons emitted from the chromophores is 
characterized as Poissonian (5^j2^,|2^], and deviation of 
the photon counting statistics from the Poisson case in- 
dicates the characteristic of fluctuations. For example, 
a super-Poissonian, photon bunching phenomenon has 
been observed in many systems with various physical ori- 
gins J30| , p2| ,^1^53|. Recently, photon counting statistics 
for a single molecule that undergoes the KA spectral dif- 
fusion process characterized by Markovian, rate processes 
has been considered (2^]2£|. In this section, we consider 
the photon counting statistics of QDs undergoing the flu- 
orescence intermittency characterized by the power-law 
process, and will show that unusual behavior in the pho- 
ton counts is obtained. 



A. Model for QD Fluorescence Intermittency 



We assume a two state model for the fluorescence in- 
tensity fluctuation of the QD: I(t) = I + or I(t) = 7_. In 
the experiments mentioned in the introduction the state 
+ is the on state while the state — is the off state and 
then 7_ = 0. We consider a more general model where J_ 
is not necessarily equal to zero, corresponding to a sin- 
gle emitter jumping between two different emitting states 
that has been observed in other situations 1541. The + 
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and — times are assumed to be mutually independent, 
identically distributed random variables. The PDFs of 
the ± times is ip±(t±). 

We consider an ensemble of N, independent, statis- 
tically identical QDs undergoing such a random pro- 
cess. Let P(n, t) be the probability of detecting n pho- 
ton counts in the time interval (0, t) from the sample. 
We use Mandel's semi-classical photon counting formula 



w n 

P(n,t) = —exp(-W), 
ni 



(3.1) 



where W is the macroscopic fluorescence intensity of the 
sample observed during the measurement time t 



W 



N 

E 

71 = 1 



N 



dt'l n (t'). 



(3.2) 



Here, w n is the contribution of the nth QD to the total 
photon counts and obtained from the fluorescence inten- 
sity I n (t) of the nth QD. £ is a coefficient which depends 
on the detection efficiency, and for simplicity, we set f = 1 
without any loss of generality. 

In the photon counting statistics of the ensemble, we 
need to consider two different averages: (i) average over 
the shot noise process d ue to the discreteness of pho- 
tons, which is denoted by (•■•) = Sr^=o(' ' ')P{ n >t): an d 
(ii) average over the stochastic process the ensemble is 
undergoing, that is, random I(t), which is denoted by 
(...) = dW(---)P(W,t), where P(W,t) is the PDF 
of the random variable W. If I n (t) is non-random and 
independent of time, I n (t) = I, the photon statistics is 
Poissonian with the mean n = W = Nit. Averaged over 
the stochastic process, Eq. ( jO| ) can be written as 

foo urn 

(P{n,t))= dWP(W,t)— exp(-W). (3.3) 



We see that {P{n, t)) is the Poisson transform of P(W, t), 
which in principle can be calculated from the statistical 
properties of the stochastic process I(t). To characterize 
the fluctuations, we use the Q parameter introduced by 
Mandel li^MM , 



(n) 



Using Eq. fl3.l| ) it is easy to show [|0|2 

n = W, 
n 2 — n = W 2 , 



(3.4) 



(3.5) 
(3.6) 



By using the fact that all QDs are statistically equivalent, 
we have 

(n) = (W) = N(w), (3.7) 
(n 5 ) - (n) = (W 2 ) = N{w 2 ) + N(N - l)(u>} 2 , (3.8) 



where we have dropped indices n in (w n ) and (w^) not- 
ing that they are independent of n. Then Q parameter 
becomes 



(w 2 ) - (w) 5 



(3.9) 



We see that Q > indicating a super-Poissonian behav- 
ior. When Q = photon counting statistics is Poisso- 
nian. 

The problem at hand is related to the Levy walk model. 
Briefly, the Levy walk model considers a test particle 
whose velocity switches randomly between two states 
v±, and the sojourn times of these two states are as- 
sumed mutually independent, identically distributed ran- 
dom variables. The PDF of sojourn times is assumed to 
decay as a power-law. In our context we may identify I± 
with the velocities of the Levy walker, and w is its coordi- 
nate. We note that there are a few variants of the Levy 
walk model (i.e., jump model, velocity model, and two 
state model) fig ]. Our model maps onto the two state 
model considered first by Masoliver et al. |35) . There are 
two technical differences between the problem at hand 
and the previous work: in our case the random walk is 
biased and asymmetric. 



B. Calculation of Q Parameter 

As shown in the previous subsection, the calculation of 
Q parameter is reduced to the calculation of the fluctua- 
tion of the random variable w since the all QDs are sta- 
tistically equivalent. We now consider the PDF and the 
characteristic function of w, P(w,t) and P(k,t), which 
are related to each other by the Fourier transform, 



P(k,t) = (exp (ik J dtl(t) 

dw exp(ikw)P(w, t). 



(3.10) 



with P(w < 0,t) = 0. We note that the characteris- 
tic function P{k,t) is mathematically equivalent to the 
relaxation function $(£) considered in the previous sec- 
tion, and therefore, can be calculated in the same way as 
before. One technical difference is that in the lineshape 
problem we have assured that the stationarity be guaran- 
teed by using the equilibrium initial condition, Eq. ( |2.7| ), 
that involves different forms for the sojourn time PDFs 
corresponding to the first transition events. Here, we use 
the same sojourn time PDFs for all transition events, 
not distinguishing between the first steps and the others. 
Thus, the process is non-stationary. 

The characteristic function is w ritte n as a sum of four 
terms in the same way as in Eq. ( ^.19| ), 



,(M) 



(3.11) 



=±n=± 
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where p m is the probability that the process begin s 
from the state to, and p + + p- = 1. In Eq. ( 3.11 ), 
Pmn(k,i) = (e lkw ) mn is the conditional characteristic 
function which is obtained by an average over paths re- 
stricted to the state to at the initial time and the state 
n at final observation time t. It can be calculated via the 
same route taken for the calculation of the conditional 
relaxation function, 



(w 2 ) - (w) 2 ~ (1 - a)p + p_ (I, 



I) 2 t 2 . 



(3.20) 



P++(k,t) = C 




i-^+(*+)iM«-) 



s_ l-'0 + (s+)r/;_(s_) 



l-^+0»+)$_(O 





-(«-) 






•5 + 


1 — -0 


+ ( S+ )^_( S _)" 



(3.12) 
(3.13) 
(3.14) 
(3.15) 



where s± = s — ikl±, and Cr x is the inverse Laplace 
transform from s to t. The equivalence between the 
charact eristic funct ion of the random photon counts w 
in Eq. ( 3.12 )-( 3.15 ) and the rel axati o n fun ction of the 
Kubo- Anderson oscillator in Eq. (2.15)-(2.18) can be seen 
explicitly if the correspondence between I± and u>± is rec- 
ognized bet ween t w o mo dels, I± <-> u>±, and f± — ip± are 
used in Eq. (|1|)-(|1|). 

Now, we consider the case where ip+ (t) and ip- (t) de- 
cay for long times as £~( 1+Q ) with < a < 1. The 
on-off intermittency of QDs corresponds to a = 1/2 if 
m om = TO ff = 3/2. In what follows, we use the Taube- 
rian theorem of the Laplace transform jllj to find the 
long time behavior of (w) and the fluctuation (w 2 ) — (w) 2 . 
We use the Laplace transforms of the sojourn time PDFs 
in the long time limit, 



$ ± (s) = 1 - A±s a 







(3.16) 



which means ip±(t) oc £~( 1+Q! ) for long times. Note 
that these sojourn time PDFs have diverging first and 
second moments. The initial conditions we choose are: 
P± = A±/(A-+A + ). 

To calculate (w) and (w 2 ), we use 



(w) = 



d(exp(ifcw)) 



(w 2 ) 



dk 

d 2 (exp(ikw)) 



dk 2 



k=0 



k=0 



Using Eq. (|3.12|)-(|3.15D we find (for < a < 1) 

(w) - {p+I++p-I-)t. 
The fluctuations are given by 



(3.17) 
(3.18) 

(3.19) 



We note that with in the context of Levy walks the be- 
havior in Eq. (|3.20| ) is called ballistic transport since the 
fluctuation (w 2 ) — (w) 2 exhibits ballistic behavior (oc t 2 ) 
instead of the normal Gaussian, diffusive behavior (oc t). 

As a second example, we consider two equal sojourn 
time PDFs which have finite first moments, (t on ) = 
(t s) — t, but diverging second moments, 

$+(s) =i>-(s) = l-TS + As a + ■■■, s->0 (3.21) 

with 1 < a < 2 and p + = p_ = 1/2. In this case we find 
in the long time limit, 



(w) 



(3.22) 



for the fluctuations we find super-diffusive behavior, 



(w 2 ) - (w) 



, A(a-l)(I + -I.y t3 _ a 



2rr(4-a) 



(3.23) 



where r(x) is the Gamma function. When a — ► 2 the 
fluctuations tend to become linear in time (for a = 2, 
Levy walks exhibit logarithmic corrections to the diffu- 
sive behavior). One can show that if a > 2, namely, the 
case when the first two moments of the ± times are finite, 
the fluctuation grows linearly with time. 



C. Discussion 



The photon statistics we find exhibits behavior which 
is very different than standard photon counting statistics. 
Specifically, Mandel's Q parameter may increase wit h 
measurement time even for long times. Using Eqs. ( [3.9] ), 
fl3~19l ), fl3~20| ), ( |3~22| ), and ( g2§ , we have 



O x { t 2 



< a < 1 

1 < a < 2 

2 < a 



(3.24) 



This is in contrast to ordinary theories of photon count- 
i ng st atistics which predict Q — > t° |5(],^6|]. From Eq. 
Q3.24 ) we see that the fluctuations are extremely large if 
compared with the standard case corresponding to a > 2. 
The mean of photon counts always increases as (n) ~ t. 
Thus, it is Q not (n) that yields insight into the under- 
lying "strange" kinetics. 

The time-independent (Q oc t°) ordinary statistics is 
a consequence of the Gaussian central limit theorem. If 
the variances of on and off time distributions are finite 
(a > 2), we expect that P(w,t) will approach a Gaus- 
sian in the long time limit, hence, the photon counting 
statistics is ordinary (i.e., p(n,t) is a Poisson transform 
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of a Gaussian, which means that photon statistics is es- 
sentially Poissonian) . On the other hand, if the variances 
of on and off time distributions diverge, standard Gaus- 
sian behavior of P(w, t) is not found even in the long time 
limit. Instead, as shown in Ref. |l6| P(w, t) will approach 
a Levy stable law (with cut-offs and 1 < a < 2). Hence, 
P(w, t) becomes very wide when a < 2 and large fluc- 
tuations occur in the photon counts. This unusual fluc- 
tuation behavior in the photon counting statistics could 
be observed in the ensemble experiments as a signature 
of the power-law blinking kinetics in QDs or other chro- 
mophores. 



IV. CONCLUDING REMARKS 



Motivated by the power-law statistics in the spectral 
diffusion and blinking kinetics in QDs observed by the 
single QD spectroscopy, we have considered two related 
phenomena characterized by power-law stochastic pro- 
cesses: lineshape and photon counting statistics. By us- 
ing the Kubo- Anderson oscillator model, we have gen- 
eralized the stochastic lineshape theory to arbitrary re- 
newal processes. Compared with the standard theory, 
we have found a variety of new phenomena in the line- 
shapes. Examples include new peaks not encountered in 
the usual KA model, and power-law narrowing behavior. 
The issue of the stationarity has been considered and the 
strong sensitivity of the lineshape to the first event in the 
stochastic trajectory was found, and it has been shown 
to be much more sensitive in the Levy sojourn time PDF 
cases with diverging first moments (0 < a < I) than the 
cases with finite first moments. 

In the second problem, we have considered the photon 
counting statistics of QDs undergoing fluorescence inter- 
mittency characterized by the power-law process, and 
showed that it exhibits unusually large fluctuations of 
photon number counts not encountered previously. The 
blinking kinetics has been mapped onto a Levy walk pro- 
cess, and due to the long time tail in the sojourn time 
PDFs Mandel's Q parameter has been shown to increase 
in time even in the long time limit. This unusual fluc- 
tuation behavior in the photon counting statistics can 
be observed in the ensemble experiments as a signature 
of the power-law blinking kinetics in QDs or other chro- 
mophores. The time-dependent behavior of Q will yield 
information on the stochastic processes the individual 
QDs are undergoing. 

Many dynamical processes in condensed phases usually 
hidden under ensemble averaging are now being revealed 
due to recent advances in single molecule spectroscopy 
p7| |6f| , and some of them turn out to exhibit "strange 
kinetics" behavior, the subject to which this issue is de- 
voted. Considering that many important systems under 
spectroscopic investigation show interesting, anomalous 
kinetic behavior, the present work or its variation will 
provide a useful theoretical framework to characterize 



"strange" kinetic behavior of complex molecular systems 
in a "usual" way. 
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APPENDIX: PROBABILITY DENSITY 
FUNCTION FOR THE FIRST SOJOURN TIME 

Although a deri vati on of the PDF for the first sojourn 
time given in Eq. (2.7) can be found in many places such 
as Refs. p2]-|45|, we present its derivation here for the 
sake of completeness. We only derive the expression of 
ijj+(t+) for simplicity. When calculating the PDF for the 
first transition event one can imagine large number of in- 
dependent stochastic processes that have been going on 
for long times before t = and have + state at t = 0. In 
general, it is not known when each of ensemble has en- 
tered the + state before t = 0, and we will call this time 
as —to- Now we define -0+(i|£o), the conditional proba- 
bility density of the transition H > — occurring at time 

t knowing that a time to has already elapsed since the 
last transition event to the + state has been made. Then 
by making use of the definition of conditional probability 
we can write as 



^+(*l*o) = 



*+(*o) 



(Al) 



In order to calculate the first transition event PDF /+(£), 
we have to average ip + (t\to) over the density ^(io), 



/+(*) 



/ °°(tt 7M*|*o)tt+(*o) 
/ °°d*o*+(*o) • 



(A2) 



With a simple change of variable it can be easily shown 
that 



/+(*) = — f°drV> + (T), 

T+ Jt 

which is Eq. (|J). 



(A3) 
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FIG. 1. A realization of the intensity fluctuation of a single QD obtained by using a two-state model. Intensity of the QD, 
I(t), fluctuates between two states, J on = 1 and J ff = 0. Statistics of on and off times are distributed with power-laws, 

fon(off )(ion(off)) ~ ^on^off j""' ' ano - m ° n = m ° s ~ were used to generate the trajectory. The lower panel is a zoom-in of the 
upper panel for < t < 2100. Note that two figures look similar on a different time scale, and in both figures long off times on 
the order of the total measurement time are observed (6600 < t < 11700 in the upper panel and 600 < t < 1370 in the lower 
panel) . 




FIG. 2. The lineshapes calculated from ^3/2^) for the stationary, (J(w), solid line) and the nonstationary (/(w), dashed 
line) cases are compared with the KA case (dot-dashed line). The correlation time is chosen as lu t = 20, 4, 0.4, and 0.08 in 
(a)-(d), respectively. The resonance frequency is chosen as ujq = 1. 
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FIG. 3. The lineshapes are calculated with the Levy PDF cases with two different values of the Levy index, a and the 
cut-off time, t c . a — 0.3 case : (a) I(lo) (stationary) and (b) I(ui) (nonstationary). a — 0.8 case : (c) I(lo) (stationary) (d) 
(nonstationary). The other parameters are chosen as lo = 1 and r± = 1. 
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FIG. 4. The power-law narrowing behavior for the lineshape with Levy PDF is shown as the Levy index a is changed, (a) 
I(uj) (stationary case) and (b) I(lo) (nonstationary case). The other parameters are chosen as r± = 0.01 and t c = 10 4 . 
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